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Abstract. The high mass X-ray binary pulsar 4U 1538-52 was observed between 
July 31 and August 7, 2003. Using these observations, we determined new orbital 
epochs for both circular and elliptical orbit models. The orbital epochs for both 
orbit solutions agreed with each other and yielded an orbital period derivative 
P/P = (0.4 ± 1.8) X 10“® yr~^. This value is consistent with the earlier measure¬ 
ment of P/P = (2.9 ± 2.1) X 10“® yr~^ at the la level and gives only an upper 
limit to the orbital period decay. Our determination of the pulse frequency showed 
that the source spun up at an average rate of 2.76 x 10“^"^ Hz sec“^ between 1991 
and 2003. 
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1. Introduction 

The source 4U 1538-52 was discovered using the Uhuru satellite (Giacconi et al. 1974). X- 
ray pulsations with a period of 529 seconds were detected from two independent satellite 
experiments: Ariel 5 (Davidson 1977) and OSO-8 (Becker et al., 1977). OSO-8 observa¬ 
tions also revealed a clear orbital modulation at 3.7 days and evidence of an eclipse lasting 
~0.51 days. The optical companion of 4U 1538-52 was identified to be the BO 1 giant QV 



Nor (Parkes, Murdin & Mason 1978). BATSE observations of this source permitted the 
construction of long-term pulse frequency and intensity histories (Rubin, Finger, Scott 
et al. 1997). In the pulse frequency history, Rubin et al. (1997) found short-term pulse 
frequency changes of either sign, and a power density spectrum of fluctuations of the 
pulse frequency derivative that is consistent with white torque noise on timescales from 
16 to 1600 days. 

From RXTE observations, Clark (2000) obtained new orbital parameters of the 
source which provided marginal evidence of orbital decay, i.e., they found Porb/Porb = 
(—2.9±2.1) X 10“® yr“^. In this work, we present new orbital epoch and pulse frequency 
measurements based on our analysis of archival RXTE observations of 4U 1538-52. 


2. Observations 

The observations of 4U 1538-52 took place between July 31 and August 7, 2003 (MJD 
52851 - 52858) and accumulated a total nominal exposure of ~ 75 ksec. The results pre¬ 
sented here are based on data collected with the Proportional Counter Array (PCA, 
Jahoda et al., 1996). The PCA instrument consists of an array of five collimated 
xenon/methane multianode proportional counters. The total effective area is approxi¬ 
mately 6250 cm^ and the field of view is ~ 1°FWHM. The nominal energy range of the 
PCA extends from 2 to 60 keV. 


3. Determination of Orbital Epoch and Pulse Frequency 

Background light curves and X—ray spectra were generated by using background estima¬ 
tor models based on the rate of very large events (VLE), detector activation, and cosmic 
X—ray background. The background light curves were subtracted from the source light 
curve obtained from the binned Good Xenon data. For 4U 1538-52, X-ray emission 
coming from the galactic ridge is only a few percent of the total X-ray emis¬ 
sion which should not affect our timing analysis (see Makishima et al. 1987), 
so we did not include an estimation of the galactic ridge emission in our 
analysis. In Figure 1, we present the background subtracted light curve. Although the 
number of active PCUs varied from 1 to 4 during the observations. Figure 1 shows count 
rates adjusted as if 5 PCUs had been active using the ’’correctlc” tool in HEASOFT 
6. For the timing analysis, we corrected the light curve to the barycenter of the solar 
system. We also corrected this barycentered light curve for binary orbital motion using 
both circular and elliptical orbital models given in Table 1 (see also Clark (2000)). Then 
a template pulse profile was created by folding the light curve into one master profile. 
Pulse profiles were also made from each of the 12 independent RXTE observations seen 
in Figure 1. 




Fig. 1 - 2-30 keV RXTE-PCA light curve of 4U 1538-52 between July 31 and August 7, 2003. 
Two 26s binned ~ 3ksec samples of this light curve corresponding to single RXTE orbits are 
presented on the upper left and the upper right. 


In order to find the pulse frequency and a new orbital epoch, we obtained 12 pulse 
arrival times through the ~ 2 binary orbits using a cross-correlation technique. In the 
pulse timing analysis, we used the method of harmonic representation of pulse profiles, 
which was proposed by Deeter & Boynton (1985). In this method, the pulse profiles for 
each orbit and the master profile are expressed in terms of harmonic series. We used 
10-term unweighted harmonic series to cross-correlate the template pulse profile with the 
pulse profiles for each RXTE observation. The maximum value of the cross-correlation 
is analytically well-defined and does not depend on the phase binning of the pulses. 

The source 4U 1538-52 has a variable pulse profile which affects the pulse timing. In 
order to estimate the errors in the arrival times, the light curve of each RXTE observation 
was divided into approximately 4-5 equal subsets and new arrival times were estimated. 
The standard deviation of the arrival times obtained from each subset of the observation 
was taken to be the uncertainty in the arrival time for that observation. 

Arrival time delays may arise from the change of the pulse frequency during the 
observation (or intrinsic pulse frequency derivative) and from the differences between 
the assumed and actual orbital and pulse parameters (Deeter, Boynton and 
Pravdo 1981), 


„ . si., ^2 c 27r(5T,r/2 aSZTO 

ocp = (po + 0V(t — to) + -zvyt — to) — OV-- - coslr, 

^ ^orbit C 


( 1 ) 


c 







where 64> is the pulse phase offset deduced from the pulse timing analysis, to is the mid¬ 
time of the observation, (j)o is the residual phase offset at to, i' is the pulse frequency 
at time to, v is the pulse frequency derivative of the source, is the orbital epoch 
when the mean orbital longitude is equal to 90 degrees, PorUt is the orbital period and 
In = 27r(t„ — Tt^/ 2 )/Porbit + 7r/2 is the mean orbital longitude at Corrected values 
of orbital and pulse parameters 5v, and v were estimated from the fits of above 

expression to pulse phase residuals. 

Table 1 presents the result of fits for both orbital and pulse parameters. 

We also reestimated the orbital epochs by varying the projected orbital radius {axlc)sini 
in the range of its uncertainty (ilcr) and found that the resulting orbital epochs are 
consistent with the best fit value at the Icr level. The error in the orbital epoch due to 
an error in {ax/c) sini may also be expressed (Deeter, Boynton and Pravdo 1981) as 


PorbO'(aa,/csini) „ nn'roj 
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where cFa,^/csini is the uncertainty in the projected orbital radius. This value is small 
relative to our error estimate for the orbital epoch (see Table 1). In Figure 2, we present 
the arrival time delay, the best fit elliptical orbit model, and arrival time residuals. 

As seen from Table 1, the orbital epochs for circular and elliptical orbital models 
agree with each other at the 1 cr level. In order to check our technique, we extracted 
observations of 4U 1538-52 done in 1997 (MJD 50449.93-50453.69) and estimated orbital 
epochs for those observations. The results agreed with the orbital epochs given by Clark 
( 2000 ). 

In Table 2, we present the orbital epoch measurements from different observatories 
and orbital cycle number (n). In Figure 3, we present observed minus calculated values 
of orbital epochs (T.n -/2 < Porbit > — < T^/2 —n< Porbit >>) relative to the constant 

orbital period (< Porbit >= 3.7228366 days). A quadratic fit to the epochs from all 
experiments yielded an estimate of the rate of period change Porb/Porb = (0.4±1.8) x 10“® 
yr“^. In Figure 4, we display the long-term pulse frequency history of the source. 


4. Discussion 

Before CGRO observations, 4U 1538-52 had been found to have a long-term spin down 
trend. A linear fit to pre-CGRO pulse frequency history gives vjv —8 x 
and a linear fit to CGRO and our RXTE result yields i>/i> ~ 1.45 x 10“^^s“^. Rubin 
et al (1997) constructed the power spectrum of pulse frequency derivative fluctuations. 
Their analysis showed that the pulse frequency derivative fluctuations can be explained 
on timescales from 16 to 1600 days with an average white noise strength of (7.6 ± 
1.6)X10“^^ (Hz s“^)^ Hz“^. A random walk in pulse frequency (or white noise in pulse 
frequency derivative) can be explained as a sequence of steps in pulse frequency with 
an RMS value of < > which occur at a constant rate R. Then the RMS variation 




Parameter 


Elliptical Orbit 


T ^/2 Orbital Epoch (MJD) 52855.0421±0.025 

a^^ sin i (lt-s)“ 56.6±0.7 


Circular Orbit 


52855.0441±0.025 


54.3±0.6 


0.174±0.015 


w“(deg) 

64± 9 


Porbit 

3.7228366±0.000032 

3.7228366±0.000032 

Epoch (MJD) 

52855.0585±0.025 

52855.0585±0.025 

^pulse (s) 

526.8551±0.016 

526.8535±0.013 

u (Hz S~^) 

(2.838±4.124)x 10“^® 

(2.241± 2.764) X 10"^® 

reduced 

1.44 

1.0 


Taken from Clark (2000) 










MJD since 4301 5.8 


Fig. 3 - The phase residuals of orbital epoch for 4U 1538-52. The orbital phases are estimated 
relative the constant orbital period (T.;r /2 — n < PorUt > — < T „/2 — n < PorUt >>), where n 
is the orbital cycle number (see Table 2). The rightmost point corresponds to the most recent 

RXTE observation of ID 80016. 

frequency derivatives can be obtained as < (A;>)^ > 1 / 2 = (5'/t)^/^Hz.s~^. As seen from 
Table 1, in our fits, upper limits on intrinsic pulse frequency derivatives are 7-10 times 
higher than the long-term spin up rates. If white noise in the pulse frequency derivative 
can be interpolated to a few days, then the upper limit on the change of frequency 
derivative obtained from a ~ 1 week observation should typically have a magnitude 

Table 2. Orbital epochs by pulse timing analysis 


Experiment 

Orbit Number 

Orbital Epoch (MJD) 

Reference 

OSO 8 

-1128 

43015.800T 0.1 

Becker et ah, 1977 

Tenma 

-457 

45517.660T 0.050 

Makishima et ah, 1987 

Ginga 

0 

47221.474± 0.020 

Corbet et ah, 1993 

BATSE 

370 

48600.979T 0.027 

Rubin et ah, 1997 

BATSE 

478 

49003.629± 0.022 

Rubin et ah, 1997 

BATSE 

584 

49398.855± 0.029 

Rubin et ah, 1997 

BATSE 

691 

49797.781T 0.022 

Rubin et ah, 1997 

RXTE 

866 

50450.206± 0.014 

Clark 2000 

RXTE 

1511 

52855.0421± 0.025 

present work 
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Fig. 4 - Pulse frequency history of 4U 1538-52. The rightmost point corresponds to most 
recent RXTE observation of ID 80016. 

that can be estimated from < >=< (A;>)^ggQ^^j^g x 15^ This value is 15/7 

- 15/10 times higher than the measured upper limit values. Therefore the 
measured upper limits on the intrinsic pulse frequency derivatives for 1 week 
are consistent with the values from the extrapolation of the power spectrum 
within a factor of a few. 

Previous marginal measurement of change in the orbital period, was (-2.9 ± 2.1) x 10“® 
yr“^ (Clark 2000), and our new value for the orbital period change, P/P = (0.4 ± 1.8) x 
10“®yr“^, are consistent with zero. These two measurements are consistent with each 
other in Icr level. 

In most of the X-ray binaries with accretion powered pulsars, the evolution of the 
orbital period seems to be too slow to be detectable. Yet there are still some such systems 
in which this evolution was measured and P/P were reported. These systems include Cen 
X-3 with (-1.8 ± 0.1)xl0“® yr“^ (Kelley et al. 1983; Nagase et al. 1992), Her X-1 with 
(-1.32 ± 0.16)X10"® yr-i (Deeter et al. 1991), SMC X-1 with (-3.36 ± 0.02) xlO"® yr"! 
(Levine et al. 1993), Cyg X-3 with (1.17 ± 0.44) xl0“® yr“^ (Kitamoto et al. 1995), 
4U 1700-37 with (3.3 ± 0.6) xlO"® yr"! (Rubin et al. 1996), and LMC X-4 with (-9.8 
± 0.7)xlO”"^ yr~^ (Levine et al. 2000). Change in the orbital period of Cyg X-3 was 
associated with the mass loss rate from the Wolf-Rayet companion star. For 4U 1700- 
37, the major cause of orbital period change was also thought to be mass loss from the 




companion star. For Her X-1, mass loss and mass transfer from the companion were 
proposed to be the reasons of the change in the orbital period of the system. 

On the other hand, for the high mass X-ray binary systems Cen X-3, LMC X-4 and 
SMC X-1, the major cause of change in the orbital period is likely to be tidal interactions 
(Kelley et al. 1983; Levine et al. 2000; Levine et al. 1993). For these three systems, orbital 
period decreases (i.e. derivative of the orbital period is negative). Our new measurement 
of orbital period change {P/P) gives the value of about —10“® yr“^ which is similar 
to the observed values of SMC X-1 and Cen X-3. Further observations can give further 
information about the orbital period change of this source. 
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